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We explore and develop the mathematics of the two multi-level ions. In particular, 
we describe some new features of quantum entanglement in two three-level trapped ions 
confined in a one-dimensional harmonic potential, allowing the instantaneous position of 
the center-of-mass motion of the ions to be explicitly time-dependent. By solving the 
exact dynamics of the system, we show how survivability of the quantum entanglement 
d ■ is determined by a specific choice of the initial state settings. 

^3 ; 

(N ■ 1 Introduction 

> 

. The ability to laser cool ions within the vicinity of the motional ground state is a key factor 

' in the realization of efficient quantum communication [Ij. Laser-cooled ions confined in 

in 
o 



an electromagnetic trap are good candidates for various quantum information processing 
^ . processes such as quantum control and quantum computing [2J. Various interactions 

models and higher order non-linear models can be implemented in this system by simply 
^ ! choosing appropriate applied laser tunings [3l H]. For more than one ion, as required in 

■ any realistic quantum logic schemes, individual addressing imposes small trap frequencies, 

whereas sideband cooling imposes high trap frequencies [5] . Experimentally, the sideband 
cooling of two ions to the ground state has been achieved in a Paul trap that operates 
in the Lamb-Dicke limit [6]. Recent advances in the dynamics of trapped ions [7j) have 
demonstrated that a macroscopic observer can effectively control dynamics as well as 
perform a complete measurement of states of microscopic quantum systems. 

On the other hand, entanglement is the quintessential property of quantum mechanics 
that sets it apart from any classical physical theory, and it is essential to quantify it in 
order to assess the performance of applications of quantum information processing |8]. 
With the reliance in the processing of quantum information on a cold trapped ion, a long- 
living entanglement in the ion- field interaction with pair cat states has been observed [S] . 
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Also, experimental preparation and measurement of the motional state of a trapped ion, 
which has been initially laser cooled to the zero-point of motion, has been reported in 
[To] . It is well known that the loss of entanglement cannot be restored by local operations 
and classical communications, which is one of the main obstacles to build up a quantum 
computer [11]. Therefore it becomes an important subject to study entanglement dynam- 
ics of qubits [m |13l [HI [151 [HI [17]. Quite recently, it has been shown that entanglement 
of two-qubit systems can be terminated abruptly in a finite time in the presence of either 
quantum or classical noises [121 [IT]. This non-smooth finite-time decay has been referred 
as entanglement sudden death (ESD) [121 [181 [T9j . At least two experimental verifications 
of ESD have been reported so far [20l [21] . Interestingly, as the entanglement can suddenly 
disappear, it can be suddenly generated, a process called entanglement sudden birth [22]. 

These recent results together with the fundamental interest of the subject, are the 
motivation for the present work: we study here the entanglement evolution of two three- 
level trapped ions interacting with a laser field taking into account the time-dependent 
modulated function. We focus our attention on the entanglement decay and generation 
due to the presence of a time- dependent modulated function in the two interacting ions and 
intrinsic decoherence. We show that entanglement of two ions in the Lamb-Dicke regime 
presents some novel features with respect to single-trapped ion. In particular, we show 
how sudden feature in the entanglement (sudden birth and sudden death of entanglement) 
occurs for the two-ion system coupled to a laser field. 

The paper is organized as follows. In section 2, we introduce the model and define 
the different parameters which will be used throughout the paper. We devote section 3 to 
an alternative approach to find an exact dynamics of the coupled system in the presence 
of the time-dependent modulated function. Section 4 is devoted to study some examples 
of entanglement evolutions which lead to different types of decay for entangled systems, 
and show how different initial state settings and time-dependent interaction affect the 
decay of the entanglement. Then we conclude the paper with some remarks on how the 
presented results could be generalized to multi-ions systems. 

2 Model 

A useful model of controlled entanglement evolution consists of two three-level ions that 
are irradiated by laser beams whose wavevectors lie in the radial plane of the trap. This 
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laser configuration excites the vibrational motion in the radial plane only {ID ionic mo- 
tion) [231 121] • Therefore, the physical system on which we focus is two three-level har- 
monically trapped ions with their center-of-mass motion quantized. We denote by a and 
the annihilation and creation operators and v is the vibrational frequency related to 
the center-of-mass harmonic motion along the direction x. Without the rotating wave 
approximation, the trapped ions Hamiltonian for the system of interest may be written 
as 

H = Hem + Hion + Hint, (1) 



where 



Hem = hva^d, 

2 

n 

i=l j=a,b,c 

,\ C(2) 



+hQ{x,t)Sg^ + hQ*{x,t)Si'J + hQ{x,t)S(^^ + hQ*{x,t)Sil\ (2) 

We denote by S^^ the atomic flip operator for the \m)i — > transition between the two 
electronic states, where S^^ = \l)ii{m\, {l,m = a,b,c). Suppose the ions are irradiated by 
a laser field of the form 53(r, t) = h~^e{i\d.p\j) exp[—i{k'r — Qt)], where e is the amplitude 
of the laser field with frequency Q and polarization vector p. The transition in the 
three-level ions is characterized by the dipole moment d and k is the wave vector of the 
laser field. Therefore if we express the center of mass position in terms of the creation 
and annihilation operators of the one-dimensional trap namely x = Ax(a^ + a), where 
Ax = {h/2vmY^'^ = r]/k is the widths of the dimensional potential ground states, in the 
X direction [rj is called Lamb-Dicke parameter describing the localization of the spatial 
extension of the center-of-mass), and m is the mass of the ions. 

By using a special form of Baker-Hausdorff theorem [3], the operator exp[z?7(a^ + a)] 
may be written as a product of operators i.e. exp{ir]{a'^+a)) = exp ^^['^^ '^]) ^^P {'^V^^) ^xp {ir]a) 
The physical processes implied by the various terms of the operator 

.M'.(a^-^a))...J=f)±^^±^^, (3, 

^ ^ n=0 m=0 

may be divided into three categories (i) the terms for n > m correspond to an increase in 
energy linked with the motional state of center of mass of the ion by {n — m) quanta, (ii) 
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the terms with n < m represent destruction of (m — n) quanta of energy thus reducing 
the amount of energy hnked with the center of mass motion and (iii) (n = m), represents 
the diagonal contributions. When we take Lamb-Dicke hmit and apply the rotating wave 
approximation discarding the rapidly oscillating terms, the interaction Hamiltonian ([2]) 
takes a simple form 

2 

Hint = ^ (nAP(t)^(aM'5i?a^ + /iAP*(t)^*(a^a)5$?a) 

i=l 

+ (hX^\t)S{a^a)S'i'^a^ + hX^^*{t)S*{a^a)SSaj , (4) 

where X^j\t) is a new coupling parameter adjusted to be time dependent. The other 
contributions which rapidly oscillate with frequency u have been disregarded. Note that 
in the Lamb-Dicke regime only processes with first and second order of t] are considered, 
while in the general case, the nonlinear coupling function is derived by expanding the 
operator-valued mode function as 

£,(a^a) = -- exp ( V ^^^l^^^at^a". (5) 

Since E^J^a^d) depends only on the quantum number a^a, in the basis of its eigenstates, 
a^a\n) = n\n), {n = 0, 1, 2, ...), these operators are diagonal, with their diagonal elements 
{n\Sk{a'' a)\n) is given by S\^\n) = —0.5e{n + fc)!)~^ra!L^(?7^) exp (— r7^/2) where L^{rj^) 
are the associated Laguerre polynomials. 



3 Solution 

3.1 Master equation 

To study the time evolution of the quantum entanglement, we will start by finding the 
general solution of the two- ion system with the Hamiltonian (Hj). The intrinsic decoher- 
ence approach [25], is based on the assumption that on sufficiently short time scales the 
system does not evolve continuously under unitary evolution but rather in a stochastic 
sequence of identical unitary transformations. Under the Markov approximation [25], the 
master equation governing the time evolution for the two-ion system under the intrinsic 
decoherence is given by 



where 7 is the intrinsic decoherence parameter. The first term on the right-hand side 
of equation (j6]) generates a coherent unitary evolution of the density matrix, while the 
second term represents the decoherence effect on the system and generates an incoherent 
dynamics of the coupled system. In order to obtain an exact solution for the density 
operator p{t) of the master equation ([6]), three auxiliary superoperators J,S and L can 
be introduced as (/i = 1) 

(7^) 



exp(jr)p(t) = Y.^H^Pm\ 

k=0 

exp (St^ p{t) = exp (^—iHr^ p(t) exp {iHt^ , (7) 
exp(Lr)p(t) = exp(-^^2)p(t)exp(-^i^2^ 



where Jp{t) = ^Hp(t)H, Sp{t) = -i[H,p{t)] and Lp{t) = -{^ /2){H'^p(t) + p{t)H^). 
Then, it is straightforward to write down the formal solution of the master equation ([6]) 
as follows 

pit) = exp( Jt) exp(S't) exp(Lt)p(0) 

00 

= 5^Mfc(t)p(0)M;(t). (8) 

k=0 



We denote by p(0) the density operator of the initial state of the system and M; 



k 



k/2 



J for all t. 

Eq. ([8]) can also be written as 



00 



H'^ exp y—iHtj exp y—'^H'^j , where the so-called Kraus operators satisfy ^ Mk{t)Ml{t) 



Pit) = e-^^* exp {-^H') {e(^«*)p(0)} exp e'^^ (9) 

where, we define the superoperator 5'mp(0) = -yHpiO)H. The reduced density matrix of 
the ions is given by taking the partial trace over the field system, i.e. p"(t) = trppit) = 
Pj^yj!j(t)|zj)(/A;|. We use the notation = \ii)®\j2), ihj = a, &, c), where |ai(2)), |&i(2)) and 
|ci(2)) are the basis states of the first (second) ion and Pij^ikit) = {ij\p°'it)\lk) corresponds 
the diagonal (ij = Ik) and off-diagonal (ij 7^ Ik) elements of the density matrix p"(t). 

It is worth referring here to the fact that, Eq. iQ represents the general solution of 
the master equation ([6]) that can be used to study the effect of the intrinsic decoherence 
when the modulated function X^j\t) is taken to be time- independent. However, to study 
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the time-dependent case as seen below, it will be much more convenient to use the exact 
wave function. 



3.2 Wave function 

At any time the state vector {"^{t)) in the interaction picture can be written as 

9 oo 

\m) = Y.Y.Mn,t)\^,). (10) 

i=l n=0 

where \ipi) = \n; 0102) , 1^2) = \n + 1; 0162) , l^/'s) = \n + 1; 0102) , IV'4) = \n + 1; 6102) , l^p^) 
\n + 2]bib2), iV^e) = \n + 2;biC2), l^y) = |n + 1; 0102) , \ips) = \n + 2]Cib2), \ipg) = 
|n + 2; C1C2) , where \n) corresponds to the initial state of the vibrational motion. We 
denote by Ai{n, t) the probability amplitude for the ion to be in the state In general, 
the dynamics of the probability amplitudes Ai{n,t) is given by the Schrodinger equation, 

^^^ = 2^C',,AKn,t), (11) 
1=1 

where Cjk = \ Hint \4'k) ■ These equations are exact for any two three-level system. 

In order to consider the most general case, we solve equation (fTTj) . by assuming a new 

9 

variable G{n,t) as G{n,t) = ^ XjAj(ri, t), where Xi = 1, which means that 

i=l 

^^^^ = 7(t)/5l (^A,{n,t) + ^^P'Mn,t)^ , (12) 

where X^j\t) = "y{t)Xj and /5j are given by (3i = X2C21 + X3C31 + X4C41 + XjCji, and (3[ = 
Pi/Pi, where P2 = Cii+a;5C5i+a;8C8i, = Cis + xqCqs+xqCqs, (3^ = C14 + X5C54+X6C64, 

/^s = 2:26*25 + a;4C'45, /5g = XsCsQ + X4C46, = Ci7 + XsCs7 + 3^9^*97; f^s ~ ^2C*28; 

/^g = X^C^g + XiCjg, 

Let us emphasize that in addition to the general form of equation f[T^ ). the present 
method is also suitable for any initial conditions. It is instructive to examine the formation 
of a general solution of the two three-level systems. Therefore, we use equation ffT2l) and 
seek G{n,t) such that iG{n,t) = z({t)G{n,t). This holds if = z and /3j = Xj. Some 
simple algebra gives rise to an equation that has nine eigenvalues such that the Zi [2S| are 
determined. Correspondingly, there are also nine corresponding eigenfunctions 

Gj{n, t) = Gj{0) exp (^-izj C{T)dT^ , (13) 
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where Gj{n, t) = J2 MjkAkiji, t), and Mji = p^e^i + J2 ^J^xi- We denote by Cxi mutually 

k=l i=l 

orthogonal unit vectors, given by e^^ = {6ii, 6i2, ^is, 5i4, ^is, ^ie, Sir, ^is, Sig), Sij = 1 if i = j 
and 6ij = if i 7^ j and p = (1,1,1,1,1,1,1,1,1), Xi = {xii,Xi2,Xi3,Xii, Xi5,XiQ, x„, 
Xis,Xig). Then, one can express the unperturbed state amplitude Ai{n,t) in terms of the 
dressed state amplitude Gj{n,t) as follows 

A{n,t) =Y,M-^G,{n,t) = ^M^iG,(0)exp (~zz, faT)dT) . (14) 
j=i j=i \ Jo J 

We have thus completely determined the exact solution of a two three-level system in the 
presence of the time-dependent modulated function. 



4 Entanglement evolution 

Good measures of entanglement are invariant under local unitary operations. They are 
also required to be entanglement monotones, that is, they must be non-increasing under 
local quantum operations combined with classical communication [271 l2Hl [29|, |30l [3T|, |32] . 
computable entanglement measures for a generic multiple level system do not exist. For 
two-qubit systems, concurrence can be used to compute entanglement for both pure and 
mixed states [33]. Rungta et al [M] defined the so-called /-concurrence in terms of a 
universal-inverter which is a generalization to higher dimensions of two spin flip operation, 
therefore, the concurrence in arbitrary dimensions takes the form 

(t) = ^/{^P\SN,0SN,m{^P\m. (15) 

Another generalization is proposed by Audenaert et al [53] by defining a concurrence 
vector in terms of a specific set of antilinear operators. As a complete characterization of 
entanglement of a bipartite state in arbitrary dimensions may require a quantity which, 
even for pure states, does not reduce to single number [36], Fan et al. defined the concept 
of a concurrence hierarchy as — 1 invariants of a group of local unitary for A^— level 
systems [37] . 

Now, we suppose that the initial state of the principle system is given by 

00 

1^(0)) = (cos^|ai,62) +sin^e'^|6i,a2)) ^^^^ 

n=0 
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Figure 1: The evolution of the concurrence (t) as a function of the scaled time Xit and 
6. The parameters are n = 5, A2/A1 = 0.01, and = 0. 

where 6 G [0,27r], G [0,7r] and = (n|V'(0)). 

In absence of intrinsic decoherence, the concurrence (1151) will be calculated using Eq. 
( ITOl) . while to study the effect of 7 we use equation ([9]). In figures 1 and 2 the dependence 
on 6 and the scaled time Ait, in the Lamb-Dicke regime, of the concurrence [t) is 
illustrated. In the typical experiments at NIST [39], ^Be'^ ions are stored in a RF Paul 
trap with a secular frequency along xofz//27r~ 11.2 MHz, providing a spread of the 
ground state wave function of Aa; ~ 7 nm, with a Lamb-Dicke parameter of r/ ~ 0.202. 
With these data we find e ~ 0.01, so they can be considered as small parameters. We 
see that the concurrence (t) exhibits some peaks whose amplitude decreases as the 
interaction time increases. The main consequence is that we can select a given value of the 
superposition parameter 6, or rather a specific instant time to obtain strong entanglement 
or disentanglement. 

The entanglement sudden birth phenomenon is observed in figure 1. Entanglement 
is not present at earlier times, and suddenly at some finite time an entanglement starts 
to build up. Also, from figure 1, it is clearly seen that the value of first local maximum 
significantly exceeds the second local maximum when the two ions start from a superpo- 
sition state. Recall that the entanglement attains the zero value (i.e., disentanglement) 
when the trapped ions start from either |ai,62) or |&i,a2) states, while strong entangle- 
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ment occurs when the inversion is equal to zero, i.e, the two ions start from a maximum 
entangled state, such as IV'^b(O)) = (ki; ^2) + \bi, 02)) /V^- In other words, for the initial 
entangled state, there are some intervals of the interaction time where the entanglement 
reaches its local maximum and drops to zero (entanglement sudden death). 




Figure 2: The same as figure 1 but n = 15. 

To gain more insight into the general behavior of the quantum entanglement evolution, 
we plot in figure 2 the time evolution of the generalized concurrence (t) for a larger value 
of the field intensity, n. It can be seen that the dynamics is strongly modified where a 
smooth decay of the entanglement is observed. In contrast to the dynamics with the small 
n, where the sudden death of entanglement-like feature has been observed after the first 
local maximum of the entanglement (see figures 1 and 2), here the entanglement smoothly 
decays and vanishes as the time increased further. Note that the effect of the superposition 
parameter 6 on the entanglement distributions in both figure 1 and 2 is similar and shows 
symmetry around = while the local maxima correspond to 6' = {n = 1,3,5, ...). 
Obviously from these figures, if ^ = (n = 0,1,2,...) i.e., the system starts from a 
separable state, the entanglement is always zero. On the other hand, if the system starts 
from an entangled state and because of the existence of further system parameters, the 
entanglement decay occurs with small or large values of n. Moreover, comparison of 
figures 1-2 shows that the entanglement decay takes a longer time to reach zero in the 
case where large n is considered. We hence come to understand that considering a large 
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intensity of the initial state of the field, can be used positively in preventing entanglement 
sudden death or delaying the disentanglement of the two ions. A physical explanation 
of why is the intensity of the field playing such a role in the disentangling process of 
the two ions, is that the strong field providing some sort of shielding mechanism to the 
decoherence effectively induced by the trace over the motional degree of freedom. 




Figure 3: The concurrence as a function of the scaled time and decoherence parameter. 
The other parameters are the same as figure 1. 

Quantum coherence and entanglement typically decay as the result of the influence of 
decoherence and much effort has been directed to extend the coherence time of systems 
of interest. However, it has been shown that under particular circumstances where there 
is even only a partial loss of coherence of each ion, entanglement can be suddenly and 
completely lost [12j. This has motivated us to consider the question of how decoherence 
effects the scale of entanglement in the system under consideration. Once the decoherence 
taken into account, i.e., 7 7^ 0, it is very clear that the decoherence plays a usual role 
in destroying the entanglement. In this case and for different values of the decoherence 
parameter 7, we can see from figure ([3]) that after switching on the interaction the entan- 
glement function increases to reach its maximum showing strong entanglement. However 
its value decreases after a short period of the interaction time to reach its minimum. 
The function starts to increase its value again however with lower local maximum values 
showing a strong decay as time goes on. Also, from numerical results we note that with 
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the increase of the parameter 7, a rapid decay of the entanglement (entanglement sudden 
death) is shown [15], [38] . 




Figure 4: The same as figure 1 but the modulated function ((t) = sech(t/2r) . 

The above results and connections are very intriguing, and naturally lead us to ask 
what is the role played by the modulated function in obtaining these associated phenomena 
of the entanglement. In order to answer this question, we consider the modulated function 
({t) to be time- dependent of the form ({t) = sech(t/2r) [101 SI]- In this form the coupling 
increases from a very small value at large negative times to a peak at time t = 0, to 
decrease exponentially at large times. Thus, depending on the value of r and the initial 
time to, various limits such as adiabatically or rapidly increasing (for to < Xt < 0) or 
decreasing (for < At < to) coupling can be conveniently studied. This allows us to 
investigate, analytically, the effect of transients in various different limits of the effect 
of switching the interaction on and off in the ion-field system. The vanishing of the 
interaction at large positive times leads to the leveling out of the inversion. 

It should be noted that the time dependence specified in ((t) is one of a class of 
generalized interactions that may offer analytical solutions. It is evident from figure 4 
that the entanglement sudden birth phenomenon is more clearly observed for the time- 
dependent interaction. In this case, the entanglement starts suddenly to build up at later 
times compared with the time-independent case. At this point the entanglement from zero 
evolves to its local maximum value and then oscillates with lower local maximum followed 
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by smooth decay (see figure 4). Although increasing the field intensity leads to strong 
entanglement (maximum value of entanglement), however the local maximum values of 
the entanglement also vary and occur for some short periods of the interaction time. This 
indicates that in a regime where coherent state is considered, the underlying states are 
highly entangled. Consequently, the presence of the time-dependent modulated function 
increases the number of oscillations and delays the disentanglement. All these results 
confirm the possibility of a practical observation of time-dependence of the modulated 
function effects for prolonging time for the disentanglement. Our results here show the 
important role played by the modulated function in the entanglement dynamics which is 
crucial for the onset of either entanglement sudden death or sudden birth in the trapped 
ion systems considered in this paper. In addition, our results suggest that the analytical 
results presented here, could be attained for different configurations of any two trapped 
ions systems [32]. 

The remaining task is to identify and compare the results presented above for the 
entanglement with another accepted entanglement measure such as the quantum relative 
entropy. Eisert and Plenio P3] have raised the question of the ordering of entanglement 
measures. It has been proven that all good asymptotic entanglement measures are either 
identical or fail to uniformly give consistent orderings of density matrices [S] . One of the 
best understood cases is entanglement measure defined in terms of the quantum relative 
entropy. More explicitly, for the entangled states p{t) the quantum relative entropy is 
defined by the following formula as the distance between the entangled state p{t) and 
disentangled state tr^p(t) ® trBp{t) e&{Hi® H2) US iH] 

Ip (P^,pf) = trp(t) (log p(t) - log(tr^p(t) ® tr^pm, (17) 

where pf- = trs {pit)) and pf = tr_4 {pit)) , A{B) refers to the first (second) ion. Note 
that if the entangled state p{t) is a pure state, S{p{t)) = and then S{tr_^p{t)) = 
S{trisp{t)), which means that we have Ip {pf,pf) = 2S{trisp{t)). One, possibly not very 
surprising, principal observation is that the numerical calculations corresponding to the 
same parameters, give nearly the same behavior with different scales. This means that the 
entanglement measured by either quantum relative entropy or concurrence measures gives 
rise to qualitatively the same results. We must stress, however, that no single measure 
alone is enough to quantify the entanglement in a multilevel system. 
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5 Conclusions 



In summary, we have derived an intuitive extension of the standard quantum model of two 
three-level trapped ions interacting with a laser field to include the time-dependent mod- 
ulated function and intrinsic decoherence. This study reveals that the time-dependent 
modulated function can be used for generating either entanglement sudden death or sud- 
den birth depending on a proper manipulation of the initial state setting. We note that 
the existence of entanglement sudden death reveals a fact that the non-interacting and 
non-communicating two ions can abruptly lose their entanglement. Also, it will be very 
interesting to extend these results to the case of mixed states in the presence of the de- 
coherence. We hope the presented results can be useful for the ongoing theoretical and 
experimental efforts in multi-levels particles interaction. 
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